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Abstract. Let du be a measure in R d obtained from adding a set of mass 
points to another measure dfi. Orthogonal polynomials in several variables 
associated with du can be explicitly expressed in terms of orthogonal polyno- 
mials associated with d/x, so are the reproducing kernels associated with these 
polynomials. The explicit formulas that are obtained are further specialized 
in the case of Jacobi measure on the simplex, with mass points added on the 
vertices, which are then used to study the asymptotics kernel functions for dv. 



Let dfj, be a measure on R with all finite moments and we assume that <i/i is 
positive definite in the sense that J Kd p 2 (x)d/u. > for every p G Il d , p ^ 0, where 
II d denotes the space of real polynomials in d- variables. Let (•, ■) ^ denote the inner 
product defined by 



Then orthogonal polynomials of several variables with respect to (p, q) M exist. Let 
N > 1 be a positive integer and let £i, £2, • • • , £zv be distinct points in R d . Let A be 
a positive definite matrix of size N x N. We define a new inner product (•, ■)„ by 



(1.2) (p, q) v := (p, q)„ + (p(fc),p(6), ■ ■ ■ ,p(6r))A(g(£i), g(6), • • ■ , q(M) 



where the superscript tr indicates the transpose, which can be defined via an integral 
as in (jl.ip against a measure dv that is obtained from adding N mass points to 
dfi. A typical example is when A is a diagonal matrix with positive entries. The 
purpose of this paper is to study orthogonal polynomials with respect to the new 
inner product (■,■)„■ 

In the case of one-variable, the first study of orthogonal polynomials for measures 
with mass points was carried out, as far as we know, by Uvarov ([9]), who gave a 
short discussion on the case of adding a finite set of mass points to a measure and 
showed how to express the orthogonal polynomials with respect to the new measure 
in terms of those with respect to the old one. The problem was later revitalized by 
A. M. Krall ([7]), who considered orthogonal polynomials for measures obtained by 
adding mass points at the end of the interval on which a continuous measure lives. 
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1. Introduction 



(1.1) 
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The case of Jacobi measure with additional mass points at the end of [—1, 1] was 
studied in [5] , where explicit formulas of orthogonal polynomials were constructed. 
For Jacobi weight with multiple mass points, it is possible to study asymptotic 
properties of orthogonal polynomials [4]. In the case of several variables, however, 
only the case of N = 1 has been studied [3] . 

Our main results contain explicit formulas that express orthogonal polynomials 
and reproducing kernels with respect to (■,■)„ in terms of those with respect to 
(•,•)„• These results are stated and proved in Section 2. As an example, we consider 
the case of Jacobi weight function on the simplex in M. d , with mass points added at 
the vertices, for which our formulas can be further specified and expressed in terms 
of the classical Jacobi polynomials. The result is then used to study the asymptotic 
expansion of the Christoffcl functions with respect to (•, •} . 

2. Orthogonal polynomials for measures with mass points 

We start with a short subsection on necessary definitions, for which we follow 
essentially [2], and prove our main results in the second subsection. 

2.1. Preliminary. Through this paper, we will use the standard multi-index no- 
tation. Let No denote the set of nonnegative integers. For a multi-index a = 
(a±, . . . , ad) G Nq and x = (xi, . . . , Xd) € K d , a monomial in d variables is defined 
as x a = x" 1 ■ ■ ■ x^ d . The integer |a| = a\ + ■ ■ ■ + ad is called the total degree of x a . 
We denote by V% the space of homogeneous polynomials of degree n in d-variables, 
P„ := span{x Q : |a| = n}, and denote by 11^ the space of polynomials of total 
degree at most n. The collection of all polynomials in <i-variables is H d . It is well 
known that 

„ri fn-\-d\ , , _j (n + d — 1 \ j 
dimLT^ = f J and dirnT^ = f J ■= rt 

Let (•, -) M be the inner product defined in (jl.ip . A polynomial p £ 11^ is orthog- 
onal with respect to (jl.ip if 

(p, q )^o, Vgen^. 

Our assumption that dfi is positive definite implies that orthogonal polynomials 
with respect to (•, -) M exist. Let us denote by V d the space of orthogonal polynomials 
of total degree n. It follows that dimV^ = r d . Let {Pa}\a\=n denote a basis of 
V d . It is often convenient to use vector notations introduced in [6] and [10]. Let 
{ai,a2, . . • ,a r d} be an enumeration of the set {a S Nq : |a| = n} according to a 
fixed monomial order, say the lexicographical order or the reversed lexicographical 
order. Then the basis {P™ }i a |=n can be written as 



J / pn pn pn \ 

n y <x\ ' a.2 3 ' * * 3 a r d J 



We will treat V n both as a set of functions and as a column vector of functions. As 
column vectors, the orthogonality of {P™. } can be expressed as 



Jtr \ 
" ml fi 



if n =/= m, 
if n = m, 



where the superscript denotes the transpose (so that P tr is a row vector) and H n is 
a matrix of size x r^, necessarily symmetric, and in fact a positive definite matrix 
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by our assumption on d\i. For convenience, we shall call the system {P n }^L = 
{P™ : \a\ = n, n = 0, 1, . . .} an orthogonal polynomial system (OPS). If H n is the 
identity matrix, then {P™ : \a\ = n} is an orthonormal basis for V% and the OPS 
is called an orthonormal polynomial system. 

Likewise, we can write x" := {x a : \a\ — n} = {x ai ,x a2 , . . . ,x r «} and regard 
it as a column vector. Since each element in P„ is a polynomial of degree n, it can 
be written as a sum of monomials, which, in vector notation, becomes 

n 

P„ = Gj,n xJ ? where G . n S M r d xr d, 

3=0 

in which M pxq denotes the set of real matrices of size p x q. In particular, G n ^ n 
is a square matrix and it is necessarily invertible since P„ is a basis of V„. We call 
G n:n the leading coefficient of P„. 

With respect to dfi, the reproducing kernel of V^f, denoted by P n (dfi; x , y) , is 
denned by (P n (d/z; x, OjP)^ = P{ x )i P € V,f. In terms of a basis P n of V^, it satisifes 

P n (dfi; x, y) = ¥l{x) H- 1 P„(y) with H n = (P„, P„)^. 

Similarly, the reproducing kernel of 11^, denoted by K n (d/i; x,y), is defined by 
(K n {dn\ x, ■),p) fJ/ =p(x), p £ 11^, and satisfies 

n 

K n (dn; x,y) = } j Pj{ d n; x, y), n > 0. 

3=0 

Since the definitions of P n (dfj,; x, y) and K n (dfi; x, y) are independent of the choice 
of a particular basis, (see [2j Theorem 3.5.1]), it is often more convenient to work 
with an orthonormal basis. The kernel K n (d)j,;x,y) plays an important role in 
studying Fourier orthogonal expansions, as it is the kernel function of the partial 
sum operator. The reciprocal of K n (d/j,; x, x) is called Christoffel function, denoted 
by A„(x), and it satisfies 

K(x) := -— -J = inf f \P(y)\ 2 dfi(y). 

K n (dfi;x,x) p(x)=i,Pen£ J U d 

2.2. Main results. Our goal is to study orthogonal polynomials with respect to 
the inner product (•,•)„ defined in (|1.2p . Let us recall that A is a given positive 
definite matrix of order N and {^1,^2, ■ ■ ■ ,£,n} is a set of distinct points in M. d . 
Introducing the notation 

p(O = {p(6).p(&),...,p(60}, 

and regarding it also as a column vector, we can then rewrite the inner product 
(•,•)„ in (O) as 



(1-2') M, = <P,<Z>„+P(O tr Aq(£), 

where (•,■)/* denotes the inner product defined in In the case that A is a 

diagonal matrix, A = diag{Ai, . . . , Xn}, the inner product (•, •) takes the form 

N 

(2-1) (p, q) v = (p, + ]T y/^.U^ji. 

3=1 

Our first result shows that orthogonal polynomials with respect to (p, q) v can be 
derived in terms of those with respect to (p, q)^. The statement and the proof of 
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this result relies heavily on the vector-matrix notation. To facilitate the study, we 
shall introduce several new notations. 

Throughout this section, we shall fix P„ as an orthonormal basis for associated 
with dfi. We denote by P n (£) the matrix that has P n (£i) as columns, 

(2.2) P„(0 := (P n (a)|Pn(6)l • • • |Pn(6v)) G M r , xN , 

denote by K„_i the matrix whose entries are K n -i(d[i; £j, 



(2.3) K n _x (if n _ 1 (d W 6,^))i W I . =1 G X 



N x N i 



and, finally, denote by K„_i(£,a;) the vector of functions 

(2.4) K n _i(£,a;) = £i, x), if„_i(d^; £ 2 , x), . . . , K n ^i(dn; £jv, x)} , 

which we again regard as a column vector. 

From the fact that K n (dfi; x, y) — K n ^i(dfi; x, y) = P n (dfi] x, y), we have imme- 
diately the following relations, 

(2.5) Pn(OP»(aO = K n fox)-Kn-i(t,x), 

(2.6) P»(OPn(f)=K n -K n _i, 

which will be used below. Let In denote the identity matrix of order N. 
Lemma 2.1. The matrix In + AK„_i is invertible. 

Proof. First we show that the matrix K n _i is positive definite, By the definition 
of K n (dn; •, •), for every c 6 R N , c ^ 0, we have 

N N 

C tr K n _ lC = E C*^a(6)^ate) = E |E C ^(0) 

|a|<n— 1 i,j=l |ce|<n-l j=0 

so that K n _i is positive definite. The matrix A is also positive definite, by as- 
sumption, so that it is invertible. Since A -1 (In + AK„_i) = A -1 + K n _i, we see 
that it is positive definite as well, hence invertible. Consequently, In + AK„_i is 
invertible. □ 

We are now ready to state and prove our first main result. 

Theorem 2.2. Define a polynomial system {Q n } n >a by Qa(x) := Po(x) and 

(2.7) Q n (x) = P n (z)-P n (0(/jv+AK n _ 1 )- 1 AK n _ 1 (£,x), n > 1. 

Then {Q n }n>o is a sequence of orthogonal polynomials with respect to (•, •) defined 
in (II. 2|) . Conversely, any sequence of orthogonal polynomials with respect to (jl.2[) 
can be expressed as in l|2.7[) . 



2 

>0, 



Proof. Let us assume that {Q„} n >o is an OPS with respect to (•, -) u and Q n has the 
same leading coefficient as P„, which implies, in particular, that Qo is a constant 
and Qo = Po- We show that Q n satisfies (|2.7|l . By the assumption, the components 
of Q n — P„ are elements in for n > 1. Since {P ra }ri>o is a basis of II , we 

can express these components as linear combinations of orthogonal polynomials in 
Po,Pi, . . . ,Pn-i- I n vector-matrix notation, this means that 



n-l 
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where M™ are matrices of size rt x rf. These coefficient matrices can be deter- 
mined from the orthonormality of P„ and Q n . Indeed, by the orthogonality of Q n , 
(Q„,Pj)^ = for < j < n — 1, which shows, by the definition of (-, •) and the 
fact that Fj is orthonormal, 

M? = (Q„,P}% = -Q„(O tr AP„(C), 

where P n (£) is defined as in (gj) and Q„(f) - {Qn(£i)|Qn(6)| ■ ■ • |Qn(6v)} in the 
analogous matrix with Q n (£j) as its column vectors. Consequently, we obtain 

n-i 

(2.8) Q„(a;) = P„(x) - ^ Q n (£) AP* r (£) Wj(x) 

3=0 

= P„(a;)-Q„(e)AK„_ 1 (^a;). 

where the second equation follows from the relation (|2.5p , which leads to a tele- 
scoping sum that sums up to K n _i(£, a;). Setting x = £j, we obtain 

Q„(e i )=Fnte)-Qn(e)AK n _ 1 (e,c i ), 1<*<^V, 

which leads to, by the definition of K n _i at (|2.3p . that 

Qn(0 = P„(0-Qn(£)AK n _!. 

Solving for Qn(£) from the above equation gives 

(2.9) Q„(0 - P„(0(/jv + AK^)" 1 , 

Substituting this expression into (|2.8p establishes (|2.7p . 

Conversely, if we define polynomials Q„ by l|2.7[) . the above proof shows that 
Q n is orthogonal with respect to (•, •)„. Since Q n and P n have the same leading 
coefficient, it is evident that {Q n }n>a is an OPS in H d . □ 

Let {Qn} n >o be an OPS with respect to (|1.2[) as in Theorem 12. 21 In general, Q n 
is not orthonormal. We denote, in the rest of this section, 

Then H n is a positive definite matrix. It turns out that both H n and H~ l can be 
expressed in terms of matrices that involve only {Pj}j>o- 

Proposition 2.3. For n > 0, 

(2.10) H n = I K + P n (0(lN + AKn-O-^K), 

(2.11) ff" 1 = I ri - P„(0(i"iv + AK„)- X AP£(£). 

Proo/. Since P„ is orthonormal, (P„,P£) M = I r d. From (|2~7| and ([2~9]) we obtain 

H n =<Q„,Q£>„ = (Q»,F£>„ = (Q„,n% + Q„(0AP t „ r (e) 

=/ r{ l + P»(0(JjV + AKn-i)-^^), 

which proves (|2~10|1 . In order to establish IpTTT]) . we need to verify that 

H n (I r 4 - Pn(0(lN + AK^AP^)) = Ir^ 

which, by (|2.10|) and after simplification, reduces to the following equation, 

(2.12) (I N +AK n _ 1 )- 1 A P n (O tr Pn(0(lN + AK,,)- 1 

= (I N + AVi)- 1 - (I N + AK*)" 1 . 
Using (|2.6[) . the above equation can be verified by a simple computation. □ 
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Our next result gives explicit formulas for the reproducing kernels associated 
with (•, •) , which we denote by 

n 

P 3 (dv; x, y) := Qj(x) Hj 1 Q^y) and K n (dv; x, y) := ^ P^dv; x, y). 

3=0 

Theorem 2.4. For j > 0. 

(2.13) Pjidv-x.y) = P j (d^x,y)-Kj(^x)(I N + A K,)- 1 A y) 

+ K t J_ 1 (£,x) (I N + AKj_ 1 ) _1 AKj_ 1 (^,y). 

Furthermore, for n > 0, 

(2.14) K n (dv;x,y) = K n (dii;x,y) - K«(£,x) (I N + AK„) _1 AI„(^,i/). 

Proof. Since A^ 1 (In +AKj_i) = A -1 + K 3 _i is a symmetric matrix, so is (In + 
AKj.i)- 1 A. Using this fact, it follows from (j!T7]) and ([2~TT]) that 

ffr 1 =p*(z) - ¥j(x)Pj(0(I - AK^-^P^O 

- K^, z) (/at IAKh)- 1 A P} r (0 

- K$U(£, x) (J w +AK j . 1 )~ 1 A PjiOP&Xl + AK^PJ^), 
which simplifies to, upon using (|2.12[) and (|2.5|) . 

Qj(x) HJ 1 =Pj(x) Pj(x)P,(0(I - AK.O-^P^O 
- KJ_^, x) (I N +AK h )-' A P?(£) 
=Pj( a; ) - K»(£, z) (J* + AK,)- 1 A PJ(0- 
Using again (|2.7[) and (|2.5p . we then obtain 
Qj(x)Hr 1 Q j (y)=Fj(x)F j (y) 

- [K* r (e, a:) - Kj^x)} (I N +AK j _ 1 )- 1 AK j _i(f, y) 

- K^, x) (Jjv + A K,)- 1 A [Kj(£, y) - K^fo y)] 

+ Kj^z) (Jjv + AKj)- 1 A P7(0Pj(0(^v +AK H )- 1 AK J '-i({, y), 

which simplifies to (|2. 13[) upon using the identity (|2. 12[) . 

Finally, summing over (|2. 13[) for j = 0, 1, . . . , n, we obtain (|2. 14[) . □ 

The result in this section can be extended without much difficulty to mass points 
with derivative values. To be more precise, let d a — 9" 1 ■ ■ • d^ d , where di = 
and for cti £ Nq, i = 1,2, ... ,N, define 

(2.15) D aP (0 ■= {d ai p(ii),d a - P (&), d a »p(£ N )} , 

and regard it also as a column vector. Instead of requiring £j ^ £j, we only assume 
that £j ^ £j when ai — otj. In other word, £j and £j can be the same as long as 
di 7^ oij . We then consider the inner product defined by 

(2.16) (p,q) v = ( P ,q)^ + D a p(O tr AD a q(0. 

When ai = for all i, this is the inner product in (|1.2p . Other interesting cases 
include, for example, 

N N 

(p, q) v = (p, q)^ + Y, + A i V K&) ■ V «fe)- 

3=0 3=0 
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Our results in Theorem 12.21 Proposition 12.31 and Theorem 12.41 still hold in this 
setting, but we need to replace P n in (|2.2|) . K n _i in (|2.3|) . and K n _i((,i) in (|2.4|) 

by 

p;(0 - (0 Q1 P„(£i)|3 aa P„(6)| ■ ■ ■ |d Qd P„(6v)) e A4 P * XJVj 



k;_! :=(a^ } a^ } iC_i(d/i;Ci^i))^ =1 eM NxN , 



K r*-i(£,- T ) = |9" 1 1 } iiL" n _i(d^;^i ) af),a" 1 2 j.iiL" n _i(d^;^2,af) ) • • ■ , d^K^^d/i; £ N ,x) j 

respectively, where 9?^ .f£r n (u, u) means that the derivative is taken with respect to 
u variable. 



Theorem 2.5. The results in Theorem \2.2\ and Theorem \2.J\ hold for the inner 
product defined in (|2.15|) when P n; K„_i and K„_i(£, a;) are replaced by P*, K*_ x 
and K*_ 1 (^, a;), respectively. 

The proof follows as before almost verbatim with little additional difficulty. 

3. Orthogonal polynomials on the simplex 

In this section we apply the general result in the previous section to orthogonal 
polynomials on the simplex 

T d := {x = (x u ...,x d ) € E d : Xi > 0,1- \x\i > 0} 

in M d , where |a;|i = X\ + . . . + Xd- 

3.1. Jacobi polynomials on the simplex. We consider the Jacobi weight func- 
tion 

w K (x) = x?- 1 ' 2 ■ ■ ■ x K /- 1/2 (i - ix^r^- 1 / 2 , m > 0, 

on the simplex, where 

r(H + ^i) 

is the normalization constant of W K such that w R j Td W K (x)dx = 1. Associated 
with W K , we consider the inner product on the simplex 

(3.1) (f,9)=w K j f(x)g(x)W K (x)dx, 



which plays the role of (•, •) when we deal with the settings of the previous section. 
For d = 1, W K is the classical Jacobi weight function, which has orthogonal poly- 
nomials Pn Kl ' K2 ' (2t — 1), where P^ ,h ^ is the classical Jacobi polynomial of degree 
n that is orthogonal with respect to (1 — t) a (l + t) b on [—1, 1] and normalized by 
Pn a ' b \l) — (™^ a )- We shall also denote the orthonormal Jacobi polynomials by 
pii' b \t). Evidently, p^n' b \t) — c n Pi a ' b \t), where the constant c„ is given by [HI 
(4.3.3)]. 

To state an orthonormal basis for V„ on the simplex, we follow [5J p. 47] and 
introduce the following notation. Associated with x — (x\, . . . , Xd) € R d , wc define 
by Xj the truncation of x, namely 

x = 0, Xj = (x 1 , . . . , xj), l<j<d, 



s 
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and associated with a — (ai, . . . ,a d ) € Nq and k — (ki, . . . , G we 

introduce, respectively, 

a j := (ay, . . . ,a d ), 1 < j < d, k j := {k 3 ■, . . . , K d +i), l<j<d+l. 

Then, an orthonormal basis associated with (|3.ip is given explicitly by 



where the parameters aj and 6j are given by 
and /i a is the normalizing constant given by 

(M + ^) 2H 



hi = 



nL( 2 i« j+1 i + M + 



J 2a, 



in which (a)k ■= a(a + 1) . . . (a + k — 1) denotes the shifted factorial. 

In this case, we also have a compact formula for the reproducing kernels, given in 
terms of the Gegenbauer polynomials C^, which are orthogonal with respect to the 
weight function (1 - f 2 )A-i/2 and normalized by C*(l) = ( n+2 „ A_1 ). The formula, 
first derived in Theorem 2.3], is given by 

(3.3) K n {W R \x,y) =-^— / C^fy^i yih H h y/x d+1 y d+1 t d +i) 

z ^ Jr-i.ii^+i 



[-i,i] d+1 
l 

]Jc Ki (l-t^-'dt, 



d+1 

X 
3=1 



where x d+1 = l-\x\ 1: y d+1 = l-\y\i, A := |k| + ^±1, and c Kj = $\{l-i]Yi 1 dt r . 

Let us denote the standard Euclidean basis of R d by {ei, . . . , e d }, where e, ; = 
(0, . . . , 0, 1, . . . , 0) with the single 1 in the i-th position. Furthermore, we set 
e d+ i = (0, . . . , 0) G R d . Then {ei, ei, . . . , e d+ {\ is the set of vertices of T d . 

Proposition 3.1. Let A = + For 1 < i < d + 1, we have 

(3.4, *.<*,;*,.,) - -^j-B^.^— V^-™ (2 „. „. 

In particular, we have 

M r\ ^ rw ^ 1 (A) n (A - Ki + l/2) ra t^.^,,! 

(3.5) Kn (W K ;e i ,e i ) = w - 1 — ^ + ^ , l<,<d+l. 

(3.6) ^„( W / K ;e l ,e J ) = + i ^^, l<i,j<d+l, 
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Proof. Since C^ n is an even function, for 1 < i < d + 1 we deduce from Q3.3P that 
K n (W K ;x, et) =^ J ^ CK^U) (1 - i*)"*" 1 dt t 



l 

cL(V^u) (i - to*- 1 (i + u)^ du 



2 d+i 

where is the operator defined in [2, Definition 1.5.1, p. 24] and the last 

equality comes from [21 Proposition 1.5.6, p. 27]. In particular, setting x = ej in 
(RT41) shows that 



K (W -e- e ) - — ^ p (A-K,-i/2, K ,-i/2) (11 

and, for i =/= j, 

K (W -e- e l - — ^ p(A-«*-i/2,K i -i/2)/_ 1 s 

K n {W K ,e 3 ,e % ) 2d+1 {Ki + 1/2)n ^n { i;. 

from which (031) and (32|) follow from 8, (4.1.1) and (4.1.4)]). □ 



3.2. Orthogonal polynomials on the simplex with mass points. We consider 
orthogonal polynomials on the simplex for the Jacobi measure with additional mass 
at each of the vertices of the simplex. In order to preserve symmetry, we shall limit 
ourself to the situation that every vertex has the same weight M > 0. In other 
words, we consider the inner product 

„ d+l 

(3.7) (f,g) v =w K / f{x)g{x)W K {x)dx + My^f{e i )g(e i ), M > 0. 
Jt * i=l 

In the language of the inner product (|1.2p , we assume that A is a diagonal matrix 
A = M Id+i and the inner product take the form of (|2.ip . 

We will further limit ourself to the case that k± = K% = • • • = Hd+i = ? > 0. 
Under this assumption, 

A = (d + l)(? + l/2). 

We further denote 

1 (A) B (A - q + l/2)„ 



A n := K n (W K ;ei,ei) 



2 d + l n\ (?+l/2) r 



(-1)" (A) r 



B n := K n (W K ; e J ,e l ) = — — , j ^ i. 
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As a result, we see that the matrix K„ defined in (|2.3|) is given by 

(A n B n ■ ■ ■ B n \ 

K n =(K n (W K ;e i ,e j ))f^ 1 = 



\B n B n 



■■ (A n — B n )Id+l + B n 



An) 



■ (A n — B n )Id+l + B n 



(1. 



I, 



This shows that K„ is a rank one perturbation of the identity matrix and, conse- 
quently, the inverse of the matrix Id+i + AK„ can be easily verified to be 

M 

(Id+i + AK„) _1 A 



[1 + M(A„ - B n )][l + MA n + dMBn] 



(1 + M A n +dM B n )I d+x — M B n 



As a result, we can now use Theorem 12.21 to derive an explicit orthogonal basis for 
the inner product (|3.7[) . which is given by 



,(x) =F n (x) + 



M 



d+l 



1 + M(A n -l ~ B n -1 



^P n (e i )iir n _ 1 (W«;a;,e i ) 



M 2 B n -Y 



[1 + M(A n -l - + MAn-X +dMB n -l] 

d+l d+l 

i=l i=l 

where {P ra } ra >o denotes the orthonormal polynomial system on the simplex T d given 
by (13. 2p . Furthermore, by Theorem 1 2. 4[ the reproducing kernel K n (dv\ x, y) for Ii n 
under the inner product (|3.7p is given by 

d+l 



K n (dv;x,y) =K n (W K ,;x,y) + 
(3.8) 



M 



l + M(A n -B„) ^ 
M 2 B„ 



K n(W K ; x, a) K n (W K ;y, a) 



[1 + M(A n - B n )] [1 + M A n +d M B n ] 
d+l d+l 

x Y Kn(W K ;x,e t )^2K n (W K ;y,ei). 

i=l i=l 

The explicit formula of the kernel allows us to derive a sharp estimate for the 
kernel K n (dv;x,y) from those for K n (W K ; x,y) and for the Jacobi polynomials. 
In the case of one variable (d — 1), such an estimate has been carried out in 0]. 
We shall give one result on the strong asymptotic of the Christoffcl function with 
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respect to dv on the simplex T d . For this purpose, we will need the following 
estimate of the Jacobi polynomials ([H (7.32.5) and (4.1.3)]): 

Lemma 3.2. For an arbitrary real number a and t £ [0, 1], 

(3.9) \Pt \t)\ < cn-^il - t + n -2)-(«+i/2)/2_ 

The estimate on [-1,0] follows from the fact that Pi a ' P \t) = (-l) n P& p ' a) (-*). 

Note that (|3.9p shows that |.Pn a | < cnr x l 2 uniformly inside a compact subset 
of (—1, 1). We derive the asymptotic for the difference K n {dv\ x, x) — K n (W K ; x, x). 

Theorem 3.3. For x in T d , 

(3.10) K n{dw) -K n{ W K]X , X )= 1 r(A-, + 1/2)^ + 1/2) 



2 d + l r(A) 

£ [Pi A "^ 1/2 '^ 1/2) (2^ - 1)] 2 (1 + Oin- 1 )) 



x 

i=l 



In particular, for x in the interior ofT d , 

lim [K n (dv; x, x) - K n (W K ;x, x)} = 0, 

n — >oo 

and the convergence is uniform in any compact set in the interior of T d . 
Proof. From (|3.8|) we deduce 



M d+1 

K n (di/;x,x)-K n (W K ;x,x) = - - —j _ - ^ K n (W K ; x, af 

n n ' %—\ 

M 2 Br, 



[l + M(A n -B n )][l+M 



4A n + dMB n ] X {f K ^ X ^) 



Ma 



2 d+1 



l + M(A n -B U/ 



E [^i A "^ 1/2 ^ 1/2) (2^-i)' 



M z CiB n 



where 



[1 + M (A n - B n )] [1 + M A n +dM B n ] 
x (j2P^- ^~ 1 ^(2x l -l) ] j , 



r _ 1 (A)„ 



2 d +i (? + l/2) n " 

By the Stirling formula for the Gamma function (see [TJ (6.1.39), p. 257]), we have 

^±^=n-(l + (0) 
T\n + 1) 

as n — > oo. Consequently, it is easy to see that the following limit relations hold: 
MCI 1 r(A - <r + l/2)r(<r + 1/2) 



1 + M(A n - B n ) 2 d +i T(A) 

MB r , 



(l + OCn- 1 )) 



_ = (-l)" r(A , ^ + 1 (2) (A-2,-1) A 0(n -l)) 

1 + MAn+dMBn K ' r(? + l/2) 
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Since A - 2? - 1 = (d - + 1/2) > for d > 2 and 
/d+i \ 2 d+i 

E P^- 112 '"- 1 ^ (2 ^ - 1) < (d + 1) E i^-^As-V*) (2 ari - 1) 

\i=l / i=l 

by the Cauchy-Schwarz inequality, it follows readily that 

d+1 2 

K- n (d^x,x)-Jir n (W lt ;a;,i)=c v £ [^ A "^ 1/2, '" 1/2) (2x i - 1)] (l + O^ 1 )) , 

i=l 

where c> s is the constant 



ex.. 



i r(A-? + i/2)r(? + i/2) 



" s 2 d + l r(A) 

This is (j3~T0| . If a; is in the interior of T d , then |Pi A ~'~ 1/2 ' ? ~ 1/2) (2^-1)1 < en" 1 / 2 , 
so that K n {dv; x, x) — K n (W K ; x, x) goes to zero as n — > oo. □ 

The asymptotic of the Christoffel function for lf K was studied in [TJ], where 
most of the results were for convergence in the interior of T d . Such results carry 
over to K n {dv\ x, y) by Thcorem l3.3l In one particular case, k, — 0, the convergence 
holds for all T d as given in |12, Theorem 2.3]: 

(3.11) lim -L^K n (W -x,x) =2 d -\ x£T d , < k < d, 

n—>oo ^ n ~r a j 

where T d denotes the k -dimensional face of T d , which contains elements of T d 
for which exactly d — k inequalities in T d — {x : x\ > 0, . . . , xj+i > 0} becomes 
equalities. In particular, T d (when none of the inequalities become equality) is the 
interior of T d , and 0-dimensional face Tq is the set of the vertices. Setting k = 0, 
so that ? = and X = (d + l)/2, we see that 

pjf-*-V**-V*)(2xi - 1) = P,[ d/2 - 1/2) (2x l - 1) = {-l) n Pt 1/2 ' d/2) (l -2a*), 

which is bounded by cn -1 / 2 whenever 1 — Xi > e > 0. It follows then that 



x £ T d , k> 

1 X £ Tn, 



n J t=i {± xkz j-o, 

upon using the fact that P^°'^(l) = (™^ a ). By (|3.10p . we then end up with the 
following corollary. 

Corollary 3.4. For k = 0, 

1 r . , |2 d - fe , fc>0, 

hm - , K n (dv. x, x) — < , " 

(«+«*) »v » » ; | 2 rf + P d , xe^ 

w/iere E d = c (d+1)/2 , = T(rf/2 + l)0F/(r(d + l/2)2 d+1 ). 

Comparing to (|3.11| . the result shows the impact of the additional mass points 
at the vertices. More generally, if Ki = ? > 0, then (|3.10p shows that 

K n {dv; x, x) - K n (W K ; x, x) =c Xs (l + ©(n^ 1 )) 

jrr 1/2 ). ^T d , 



ijfc/n+f-l/2 



), x £ Tjj, 1 < k< d- 1, 
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since, for x G T , xi = 1 only when x = ei. In particular, we see that 

lim 7n~+d\ \ K n{dv\ x, x) - K n (W K ; x, x)} = 0, x e T fe d , 1 < k < d, 

if d > 2<j — 1, whereas this limit is unbounded when x GTq. 
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